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I present recent results on the theory of QCD matter production in high energy heavy 
ion collisions and on the interactions of heavy quarks in such environment. The centrality 
and rapidity dependence of hadron production is evaluated in semi-classical approach. 
00 ■ The energy loss of heavy quarks in matter is computed. The heavy-to-light meson ratio 
^ ' (e.g., D/it) at moderate transverse momenta is demonstrated to be both sensitive to the 
density of color charges in the medium and infrared stable. 
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1. Foreword 

This conference celebrates the contributions made by Helmut Satz to the theory of 
super-dense matter. Helmut was among the very first who started to think about the 
r-| ! subject at the time of its infancy. His work on QCD matter over the years brought the field 
of statistical QCD to maturity. His strong support was vital in establishing experimental 
heavy ion programs at CERN and BNL. All of this has been already described at the 
conference by people who have a first-hand knowledge of the history of the field (in which 
the previous Bielefeld 1980 meeting was a major milestone). What I would like to add to 
these accounts is Helmut's influence on young physicists: he is, and has always been, an 
inspiration for people entering the field. Many of us, like myself, were brought into the 
field and encouraged by him. Working with Helmut is both enriching and enjoyable. I 
wish new generations of physicists will discover this for themselves. 

2. Statistical QCD: from small x to high T 

2.1. Quantum statistics at small x 

Small X physics is not usually considered to belong to the realm of statistical QCD. 
Nevertheless, especially at this conference, it is worthwhile to emphasize that the concepts 
of statistical approach provide a very useful perspective in small x physics as well. 

Let us begin by noting that a parton fluctuation with a given Bjorken x and transverse 
momentum k± inside a hadron with a momentum P has, in the Lab frame, a lifetime 
given by the uncertainty relation: 
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This shows that the partons at larger and smaller k±, live much longer than partons 
at small x and large k±. 

Let us denote hy (p = ^ip^ipjA^ the set of parton fields with x < Xq (with xq setting 
some arbitrary boundary between "fast" and "slow" fields), and by (p the set of the same 
fields, but with x > xq. Suppose that we want to compute an expectation value of some 
observable O. In doing so, we have to take account of the fact that partons at larger 
X > Xq are effectively "frozen"; this is done by employing the form familiar from the 
treatment of statistical systems with random, frozen impurities: 

where S is the action, p(0) describes the distribution of large x partons, and we have 
explicitly written down the Planck constant h. The meaning of (|]) is simple - the "frozen" 
fields are not a dynamical part of the system of the "fast" fields ip; rather, they act 
as impurities, or sources. This formulation is the basis for McLerran-Venugopalan model 
of hadron structure at small x; since glasses are among the physical systems with large 
relaxation time, one may also call such system a "color glass condensate" (see 0,0] and 
references therein). Renormalization group equations with respect to the changing scale 
Xo allow then to reconstruct QCD evolution: partons radiated by sources at larger x 
themselves become sources for radiation at even smaller x. 

At sufficiently small x and/or large atomic number of the nucleus, the density of partons 
will become very large and the system will thus cease to be dilute. What will it look like? 
What kind of dynamics will govern its properties? To address these questions, let us first 
note that for a system with large number of gluons the action is large, S ^ h. Such 
systems are appropriately described by using the semi-classical approximation. To go 
further, we need to establish the dependence of the action on the coupling constant. To 
do this, let us re-scale the gluon fields in the QCD Lagrangian as follows: — = gA^^^. 
In terms of new fields, G^^, = gG^^ = 9^v4^ — S^AJ^ + f"-^'^A^^Al, and the dependence of 
the action on the coupling constant is given by 

jj, G^^Gl, d'x. (3) 

Let us now consider a classical configuration of gluon fields; by definition, in such 
a configuration does not depend on the coupling, and the action is large, S ^ %. The 
number of quanta in such a configuration is then 

~ I ~ — P4^4, (4) 

where we re-wrote (|^) as a product of four-dimensional action density and the four- 
dimensional volume V4. 

The effects of non-linear interactions among the gluons become important when d^A^ ~ 
(this condition can be made explicitly gauge invariant if we derive it from the expansion 
of a correlation function of gauge- invariant gluon operators, e.g., G^). In momentum 
space, this equality corresponds to 

Ql ^ A^ ^ {G^fl- = (5) 
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Qs is the typical value of the gluon momentum below which the interactions become 
essentially non-linear. 

Consider now a nucleus A boosted to a high momentum. By uncertainty principle, the 
gluons with transverse momentum are extended in the longitudinal and proper time 
directions by ~ l/Qs; since the transverse area is vri?^, the four-volume is V4 ~ 'kR\/Q1. 
The resulting four-density from (HI) is then 



1/4 '^^A 



where at the last stage we have used the non-linearity condition (H), p4 ~ Q^. 

Identifying the number of gluons in the infinite momentum frame with the gluon struc- 
ture function Ng ~ xG{x, Ql), we arrive at the condition 

^2 ^ xGa{x,Q1) , , 



originally derived in [^n]| as the criterion for "parton saturation" (for a discussion of 
saturation in terms of the partons in the final state, see This simple derivation 

illustrates that the physics in the high-density regime can potentially be understood in 
terms of classical gluon fields. 

2.2. Classical QCD and particle production in heavy ion collisions 

The energy dependence of saturation scale Qs is determined by the x— dependence of 
the gluon structure function (see (|^)). In spite of significant uncertainties in determination 
of the gluon structure functions, the following observation is very important: the HERA 
data exhibit scaling when plotted as a function of variable 



g2 /x^^ 



(8) 



where A ~ 0.25 ^ 0.3. In saturation scenario, this scaling translates in the following x 
dependence of dimensionful scale Qs'- 

Ql{x) = Ql {xo/x)\ (9) 

Since the rapidity y and Bjorken variable are related by Inl/x = y, (^ leads to the 
dependence of the saturation scale on rapidity: 

Qlis; ±y) = Qlis; y = 0) exp(±Ay). (10) 

Let us now evaluate the rapidity and centrality dependences of hadron production in 
heavy ion collisions basing on this picture We need to evaluate the leading tree 

diagram describing emission of gluons on the classical level, see Fig. |l[ 

To do this, we introduce the unintegrated gluon distribution ipA^x, k'j) which describes 
the probability to find a gluon with a given x and transverse momentum kt inside the 
nucleus A. As follows from this definition, the unintegrated distribution is related to the 
gluon structure function by 
2 

/Pt 
dkl ipA{x,kl)] (11) 
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Figure 1. Mueller diagram of classical gluon 
emission. 

when pf > Ql, the unintegrated distribution corresponding to the bremsstrahlung radia- 
tion spectrum is 



M^,k^)r^^-,. (12) 
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In the saturation region, the unintegrated gluon distribution has only logarithmic depen- 
dence on the transverse momentum: 

^^{x,k^)^^; k^<Ql (13) 
as 

where Sa is the nuclear overlap area, determined by the atomic numbers of the colliding 
nuclei and by centrality of the collision. 

The differential cross section of gluon production in a AA collision can now be written 
down as 



da AnNc 1 /" ,,2 / 7 2\ / / 7\2\ 

= _ X ^ J ^f^t (^s VA{xi,kt) (Pa{x2,{p- k)t), (14) 

where Xi^2 = {Pt/ V^) Gxp{±ri), with r] the (pseudo)rapidity of the produced gluon; the 
running couphng ag has to be evaluated at the scale = max{k^, {p — k)f}. The rapidity 



density is then evaluated from (|1J) according to 



(15) 



dN 

dy ctaa 

where aAA is the inelastic cross section of nucleus-nucleus interaction. Evaluation of 
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Figure 2. Centrality dependence of charged 
hadron production per participant at differ- 
ent pseudorapidity rj intervals in Au — Au 
collisions at ^/s = 130 GeV, from |T^; the 
data are from [|15|. 



Figure 3. Pseudo-rapidity dependence 
of charged hadron production at different 
cuts on centrality in Au — Au collisions at 
^/s = 130 GeV, from the data are 



from [15 



Eqs. (|l4) and (15) leads to the following simple analytical formula [|T^, which exhibits the 
scaling properties of hadron multiplicity in nucleus-nucleus collisions: 



dN 
dy 



part 



So 



In 



Ql 



A2 



-My\ 



1 + A|y| 1 



Q 



£ Jl+X/2)\y\ 



, (16) 



with Ql{s) 



sq) {s/soY^"^. This formula expresses the predictions of high density 
QCD for the energy, centrality, rapidity, and atomic number dependences of hadron mul- 
tiplicities in nuclear collisions in terms of a single scaling function. Once the energy- 
independent constant c ~ 1 and Q1{sq) are determined at some energy sq, Eq. (|TBp 
contains no free parameters. 

The results for the Au — Au collisions at y/s = 130 GeV based on Eq . (p!6D are presented 
in Figs and ^. One can see that the agreement with the data is quite good. If it 
persists at higher energies, one may conclude that parton saturation indeed adequately 
describes the initial conditions created in relativistic heavy ion collisions. Since saturation 
provides good conditions for parton thermalization 0, we may expect that the final goal 
of producing the equilibrated QCD matter in the laboratory may be within reach. We 
thus have to look for the probes which can be used for its diagnostics. One, by now 
famous, probe of QCD matter is the heavy quarkonium [Q. Another probe is provided 
by high pt jets (see [0). I am now going to discuss a recent proposal, involving heavy 
quarks at high pt lOi 
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3. Heavy quark energy loss in QCD matter 

Let us begin by recalling the basic features of gluon radiation caused by propagation of a 
fast parton (quark) through QCD medium. As was pointed out in |T^, the accompanying 
radiation is determined by multiple rescattering of the radiated gluon in the medium. The 
gluon, during its formation time given again by (|l]) 

tform — J 2 ' (-^''') 
/v_L 

accumulates a typical transverse momentum 



.2 ^ ^2 tj^^ ^^g^ 



with A the mean free path and //^ the characteristic momentum transfer squared in a 
single scattering. This is the random walk pattern with an average number of scatterings 
given by the ratio t/orm/A. 

Combining (|TB|) and (|I^) we obtain 



describing the number of scattering centers which participate, coherently, in the emission 
of the gluon with a given energy u. For sufficiently large gluon energies, u > /i^A, when 
the coherent length exceeds the mean free path, Ncoh > 1. In this situation the standard 
Bethe-Heitler energy spectrum per unit length describing independent emission of gluons 
at each center gets suppressed: 



dW / dW \™ 1 _ a,Cn [J^ 



dujdz \duodz ) Ncoh vrcuA \ uj 




(20) 



Here Cr is the "color charge" of the parton projectile {Cr = Cp = ^2^^ ~ ^/^ the 
quark case we are interested in). 

In ( |20| ) we have substituted the characteristic ratio /x^/ A by the so-called gluon transport 



coefficient |j20| 
da 



q = p 



(21) 



which is proportional to the density p of the scattering centers in the medium and describes 
the typical momentum transfer in the gluon scattering off these centers. 

as 



The transport coefficient for cold nuclear matter was expressed in pO 



p[xG{x,Q% (22) 



with p ~ 0.16 fm ^ the average nuclear density and [xG(a;)] the gluon density in a nucleon. 
Taking ~ 0.5 and [xG'(x)] ^ 1 (at x < 0.1), yields 

gcoid ^ 0.01 GeV^ ~ 8 p. (23) 
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This estimate is an agreement with the result of the analysis of the gluon p±_ broadening 
from the experimental data on J /ip transverse momentum distributions which in the 
present notation yielded 

g= (9.4 ±0.7) p. (24) 



An estimate pO| for a hot medium based on perturbative treatment of gluon scattering 
in quark-gluon plasma with T ~ 250 MeV resulted in the value of the gluon transport 
coefficient of about factor twenty larger than (^31): 



ghot ^0.2 GeV^ ~ 20qcoid. (25) 

Multiplying ( pOD by the length L of the medium traversed,^ we arrive at the following 
expression for the inclusive energy distribution of gluons radiated by a quark: 

u < uJi = qL . (26) 



The fact that the medium induced radiation vanishes for uj > ui has a simple physical 
explanation, as according to ([19|) the formation time of such gluons starts to exceed the 
length of the medium: 



^ form 



X- J^ = J- = L- J— > L. 



Another important feature of medium induced radiation is the relation between the 
transverse momentum and the energy of the emitted gluon. Indeed, from (|r^) and 
(see also (ETF)) we derive 



kl ~ ^q io. (27) 

This means that the angular distribution of gluons with a given energy u is concentrated 
at a characteristic energy- (and medium-) dependent emission angle 

O^t^^aX". (28) 

u \^ J 

Gluon bremsstrahlung off a heavy quark differs from the case of a massless parton 
(produced in a process with the same hardness scale) in one respect: gluon radiation is 
suppressed at angles smaller than the ratio of the quark mass M to its energy E. Indeed, 
the distribution of soft gluons radiated by a heavy quark is given by 

as Cf du kjdkl n -M , . 

"^^-^r- {kl+uXr ^ ^ 

where the strong coupling constant should be evaluated at the scale determined by 
the denominator of (pQ]). Equating, in the small-angle approximation, k\_ with uj9 we 
conclude that the formula ( ^9]) differs from the standard bremsstrahlung spectrum 



asCpdudk"^: asCpdudO"^ 
dPo ~ = ^ 30 

71 to kj_ Tl UJ 0^ 

^For the sake of simplicity we assume here that the medium is static and uniform. 
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by the factor 

dPnQ = dPo- (l + ^j (31) 

This effect is known as the "dead cone" phenomenon. Suppression of small-angle radi- 
ation has a number of interesting implications, such as perturbative calculability of (and 
non-perturbative A/M corrections to) heavy quark fragmentation functions |22|,^, mul- 
tiplicity and energy spectra of light particles accompanying hard production of a heavy 
quark [^. 

In the present context we should compare the angular distribution of gluons induced 
by the quark propagation in the medium with the size of the dead cone. To this end, for 
the sake of a semi-quantitative estimate, we substitute the characteristic angle (^) into 
the dead cone suppression factor (^Tj) and combine it with the radiation spectrum ( pUD to 
arrive at 

Tf N dW as Cf 1 
where 

i^r^'^Q . (33) 

To see whether the finite quark mass essentially affects the medium induced gluon yield, 
we need to estimate the product iu for the maximal gluon energy ^ c<Ji to which the 
original distribution (pOD extends: 

M = g-/^ (I)''' ■ = {^y' , i^HQ - M^^. (34) 

This shows that the quark mass becomes irrelevant when the quark energy exceeds the 
characteristic value E^q which depends on the size of the medium and on its "scattering 
power" embodied into the value of the transport coefficient. 

Which regime is realized in the experiments on heavy quark production in nuclear 
collisions? Taking M = 1.5 GeV for charm quarks and using the values (|23|) and (^) we 
estimate 

^HQ = Vgcoid i^'/' M ^ 20 GeV i^-^j , (35) 

L 



^HQ = V^hot L^^' M ^ 92 GeV j , (36) 

for the cold and hot matter, respectively. We observe that for the transverse momentum 
(energy) distributions of heavy mesons the inequality E <C E^q always holds in practice, 
especially for the hot medium. We thus conclude that the pattern of medium induced 
gluon radiation appears to be qualitatively different for heavy and light quarks in the 
kinematical region of practical interest. 

The issue of in-medium quenching of inclusive particle spectra was recently addressed 



in p5|. The pi_ spectrum is given by the convolution of the transverse momentum dis- 



tribution in an elementary hadron-hadron collision, evaluated at a shifted value + 
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with the distribution -D(e) in the energy e lost by the quark to the medium- induced gluon 
radiation: 



da 



med 



da™" 



da"" 



dp\ J ^ ^ ' dp]_ ^"^"^ ' ' dp] 



-ip±) ■ Q{p±), 



(37) 



with Q{p±) the medium dependent quenching factor. The two facts, namely that in the 
essential region e -C p_L and that the vacuum cross section is a steeply falling function, 
allow one to simplify the calculation of the quenching factor Q by adopting the exponential 
approximation for the e- integral in (|37|): 



QiPA 



J de D{e) expj-^ -rj , C 



d 



In 



da'' 



dpi 



d\npj_ 

This integral results in the Mellin moment of the quark distribution, 

C 



Q{p±) 



D(u) 



exp 



-u I duj Niuj) e 
/o 



(3J 



(39) 



where N{uj) is the integrated gluon multiplicity defined according to (see for details) 



N{uj) = / duj 
The use of 



Qh{pa 



exp 



du' 

furnishes our final result: 



q 



vr V%± 9v^ "^V^' + Pi, 



1/3- 



(41) 



The first term in the exponent in (41) represents the quenching of the transverse momen- 
tum spectrum which is universal for the light and heavy quarks, (modulo the difference of 
the £, parameters determined by the p±_ distributions in the vacuum). The second term 
is specific for heavy quarks. It has a positive sign, which means that the suppression of 
the heavy hadron p\_ distributions is always smaller than that for the light hadrons. This 
is a straightforward consequence of the fact that the heavy quark mass suppresses gluon 
radiation. At very high transverse momenta, both terms vanish - this is in accord with 
the QCD factorization theorem, stating that the effects of the medium should disappear 
as p_L oo. How fast this regime is approached depends, however, on the properties of 
the medium encoded in the value of the transport coefficient q and in the medium size L. 

Constructing the ratio of the quenching functions, we estimate the heavy-to-light en- 
hancement factor as 



Qh{p. 



Ql{pa 



exp 



L 



q M2 



1/3- 



9^3 V^'+pI, 



(42) 



Basing on (^2|) we find [|T^ that hot QCD matter leads to a strong, factor of 2 3, 
medium-dependent enhancement of the heavy quark yields with respect to the yield of 
light quarks at moderately large pt > M. Experimentally, this effect should manifest itself 
as an enhancement of the heavy-to-light ratios such as D /ti. It will also be of interest to 
study the B/D ratio. 
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